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ABSTRACT 

The  Oriented  Scintillation  Spectrometer  Experiment  (OSSE)  on  NASA’s  Comp¬ 
ton  Gamma-Ray  Observatory  has  completed  numerous  observations  of  the  Galactic 
plane  and  Galactic  center  region.  A  principle  objective  of  these  observations  was 
to  measure  the  distribution  of  the  Galactic  positron  annihilation  radiation.  Most  of 
these  observations  provided  positive  detections  of  the  narrow  511  keV  annihilation 
line.  These  data  were  fitted  using  several  diffuse  distribution  models  representing 
various  populations  of  progenitor  objects.  The  only  model  investigated  which  is  not 
rejected  is  a  two-component  distribution  consisting  of  spheroidal  and  disk  compo¬ 
nents;  all  other  models  investigated  can  be  rejected  at  the  >  4.5cr  confidence  level. 
The  size  of  the  spheroidal  component  is  found  to  be  <  500  pc;  the  disk  component  of 
the  distribution  is  not  well  constrained  by  the  current  observations.  Adding  a  point 
source  of  emission  at  the  location  of  the  black  hole  candidate  IE  1740.7-2942  does  not 
significantly  improve  the  model  fits.  There  is  no  evidence  for  significant  time  variabil¬ 
ity  of  the  511  keV  line  flux  during  observations  of  the  Galactic  center  region.  When 
compared  to  historical  observations  of  the  Galactic  center  region  by  other  instruments, 
the  two-component  diffuse  model  suggests  a  limit  to  a  time- variable  component  of  the 
511  keV  line  flux  of  ~  4  x  10-4  7  cm-2  s_1. 

INTRODUCTION 

Numerous  observations  of  positron  annihilation  from  the  direction  of  the  Galactic 
center  have  been  performed  over  the  past  20  years  (see  Purcell  et  al.  1993;  Lingenfelter 
&  Ramaty  1989a  and  references  therein).  These  observations  have  utilized  various 
types  of  detectors  and  have  included  both  satellite  and  balloon-borne  instruments. 
Nearly  all  of  these  observations  have  resulted  in  positive  detections  of  a  511  keV 
line  characteristic  of  positron  annihilation.  The  line  centroid  and  width  have  been 
measured  using  high  resolution  Ge-based  detectors.  The  line  is  found  to  be  within 
~  0.5  keV  of  511  keV  and  to  have  a  full-width  at  half-maximum  (FWHM)  of  2.5  ±  0.4 
keV  (Leventhal  et  al.  1993).  The  intensity  of  the  line  has  been  found  to  be  positively 
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correlated  with  the  instrument  field- of- view,  indicating  there  is  a  diffuse  component 
to  the  observed  emission  (see  Teegarden  1994).  There  have  also  been  reports  of  time- 
variability  in  the  511  keV  line  flux  when  the  Galactic  center  region  was  observed  at 
different  times  by  the  same  instrument  (Riegler  et  al.  1981;  Leventhal  et  al.  1982; 
Leventhal  et  al.  1986),  suggesting  there  is  also  a  point  source  component  to  the 
observed  emission. 

The  diffuse  component  of  the  511  keV  line  emission  could  be  produced  by  sev¬ 
eral  different  mechanisms  (see  Lingenfelter  &  Ramaty  1989b  and  references  therein), 
including:  cosmic-ray  interactions  in  the  interstellar  medium,  gamma-ray  bursts,  pul¬ 
sars,  and  /3+-decay  products  from  radioactive  nuclei  ( e.g .,  56Co,  44Sc,  and  26A1)  pro¬ 
duced  by  supernovae,  novae,  or  Wolf-Rayet  stars.  Photon-photon  pair  production  in 
the  vicinity  of  an  accreting  black  hole  has  been  suggested  as  a  possible  source  of  time 
variable  annihilation  emission  (Ramaty  et  al.  1992). 

One  candidate  for  a  point  source  of  positrons  is  the  black  hole  candidate 
IE  1740.7-2942.  Observations  by  the  SIGMA  instrument  have  identified  IE  1740.7- 
2942  as  an  extremely  time  variable  hard  x-ray  source  (Churazov  et  al.  1993a).  During 
an  ~  13  hour  period  on  1990  October  13-14,  the  SIGMA  instrument  observed  a  strong 
flare  from  IE  1740.7-2942  which  exhibited  an  intense  (~  10-2  7  cm-2  s_1),  broad 
line  feature  (FWHM  ~  240  keV)  with  an  energy  of  ~  470  keV  (Bouchet  et  al.  1991). 
This  feature  has  been  interpreted  as  a  broadened,  red-shifted  511  keV  line  resulting 
from  positron  annihilation  occurring  in  an  accretion  disk  surrounding  IE  1740.7-2942 
(Sunyaev  et  al.  1991;  Ramaty  et  al.  1992).  Similar  outbursts  were  observed  by  SIGMA 
in  1991  October  and  1992  September  (Churazov  et  al.  1993b;  Cordier  et  al.  1993), 
though  a  simultaneous  observation  of  the  Galactic  center  region  by  the  OSSE  instru¬ 
ment  during  the  1992  September  outburst  did  not  detect  the  reported  excess  (Jung 
et  al.  1994).  There  is  no  direct  evidence,  however,  that  the  transient  line  feature  ob¬ 
served  by  SIGMA  is  related  to  the  narrow  511  keV  line  which  has  been  observed  since 
the  early  1970’s.  Further,  both  the  SIGMA  and  GRIP  instruments  have  provided 
95%  confidence  upper  limits  of  ~  4  X  10-4  7  cm-2  s_1  for  time-averaged  narrow  511 
keV  line  emission  from  IE  1740.7-2942  (Lei  et  al.  1993,  Heindl  et  al.  1993). 

The  Oriented  Scintillation  Spectrometer  Experiment  (OSSE)  on  NASA’s  Comp¬ 
ton  Gamma-Ray  Observatory  (GRO)  satellite  provides  the  unique  ability  to  perform 
high-sensitivity  observations  of  the  Galactic  plane  and  Galactic  center  region  to  mea¬ 
sure  the  distribution  of  positron  annihilation  radiation  and  to  search  for  time  vari¬ 
ability  of  the  emission.  The  OSSE  observations  of  the  Galactic  plane  and  Galactic 
center  region  collected  over  the  period  1991  July  12  -  1993  September  7  have  been 
used  to  investigate  the  distribution  of  the  diffuse  511  keV  line  emission  and  to  search 
for  possible  time  variability  of  the  emission.  These  data  are  found  to  be  best-fit  by 
a  two-component  distribution  model  consisting  of  disk  and  spheroidal  components. 
All  other  distribution  models  investigated  can  be  excluded  at  the  >  4.5cr  level.  There 
is  no  evidence  for  time-averaged  narrow  511  keV  line  emission  from  the  black  hole 
candidate  IE  1740.7-2942.  There  is  no  evidence  for  significant  time  variability  of 
the  511  keV  line  flux.  The  best-fit  distribution  model  is  also  compared  with  previous 
results  by  other  instruments  to  place  historical  limits  on  the  time  variability  of  the 
511  keV  line  emission. 


MODEL  DISTRIBUTIONS 

Several  diffuse  distribution  models  were  investigated  representing  various  pop¬ 
ulations  of  possible  progenitor  objects.  With  the  exception  of  the  CO  Distribution, 
which  is  based  on  Galactic  CO  observations  (Dame  et  al.  1987),  the  distributions  are 


based  on  geometric  models  for  the  volume  emissivities,  e,  throughout  the  Galaxy.  For 
each  geometric  model,  a  surface-brightness  model  was  generated  on  a  1°  X  1°  grid  in 
Galactic  coordinates  by  numerically  integrating  the  volume  emissivity  function  along 
each  line- of- sight.  The  integration  was  confined  to  a  cylinder  of  radius  R  =  20  kpc  and 
height  h  =  15  kpc,  centered  on  the  Galactic  center,  and  used  an  integration  step  size 
of  0.25  kpc  and  a  solar  galactocentric  radius  of  i?@  =  8.5  kpc.  The  surface-brightness 
models  were  normalized  to  the  total  flux  within  0.5  radian  of  l  =  0°,  b  =  0°.  The 
form  of  the  volume  emissivities  for  each  of  the  models  is  described  below. 

The  CO  Model,  originally  described  by  Leising  and  Clayton  (1985),  is  based  on 
the  results  of  CO  surveys  and  consists  of  a  radial  surface  emissivity  of  the  form: 

eco(p,z)  =  /'3exp  (1) 


where  p  is  the  galactocentric  radius  and  z  is  the  distance  out  of  the  Galactic  plane, 
both  in  units  of  kpc.  The  scale  height  of  120  pc  was  used  based  on  the  observed  scale 
height  of  interstellar  dust  and  gas  in  the  Galaxy  (Mihalas  &  Binney  1981). 

The  M31  Nova  Model,  originally  described  by  Leising  and  Clayton  (1985),  is 
based  on  the  observed  radial  distribution  of  novae  in  M31  and  is  represented  by: 
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The  Galactic  scale  height  of  300  pc  was  assumed  based  on  the  observed  scale  height 
of  novae  in  the  Galaxy  (Mihalas  &  Binney  1981). 

The  Hot  Plasma  Model  of  Skibo  et  al.  (1992),  which  represents  an  empirical 
model  for  the  distribution  of  hot  plasma  as  traced  by  the  6.7  keV  Fe  line,  was  also 
used.  This  model  is  given  by: 
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The  Visual  Luminosity  Model  is  based  on  a  model  for  the  total  star  density  in 
the  Galaxy  (Bahcall  and  Soneira  1980),  and  is  given  by: 

eVL  (p» z)  =  0-15 ed(p,  z)  +  0.00019  €f^\  (4) 

es  (,-K©  ) 

where  ed(p,  z )  is  the  volume  emissivity  for  the  disk  component  and  es(R )  is  the  volume 
emissivity  for  the  spheroidal  component  at  a  distance  R  from  the  Galactic  center. 
These  volume  emissivities  are  given  by: 
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The  Nova  Model  of  Higdon  and  Fowler  (1989)  also  consists  of  separate  disk  and 
spheroidal  components  and  is  given  by: 

£nm(^  z)  =  ed(p , z )  +  0.00557es(i2)  (5) 

where  ed(p,z)  and  es(R )  are  the  volume  emissivities  in  the  disk  and  spheroidal  com¬ 
ponents,  respectively.  These  volume  emissivities  are  given  by: 

ed(p,z )  =  exp(— 44.5  z2  —  0.297(p  —  i?®))  (5a) 
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The  Free  Spheroidal  Model  was  constructed  to  better  characterize  the  size  of  the 
spheroidal  component  of  the  distribution.  This  model  is  also  composed  of  spheroidal 
and  disk  components  and  is  given  by: 


eps  (p,  z)  =  Ad  x  ed(p,  z)  +  As  x  es(R) 


(6) 


where  Ad  and  As  are  free  parameters  of  the  fit  representing  the  amplitudes  of  the 
disk  and  spheroidal  components,  respectively.  The  volume  emissivities  of  the  disk 
and  spheroidal  components  for  this  model  are  given  by: 
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Note  that  equations  6a  and  6b  have  the  same  form  as  the  spheroid  and  disk 
components  of  the  Visual  Luminosity  Model  given  in  equations  4a  and  4b,  with  the 
spheroid  core  size  ( Rcore )  and  scale  size  ( Rscaie )  as  additional  parameters  in  equa¬ 
tion  6b.  For  this  distribution,  models  were  generated  on  a  grid  with  various  values 
of  the  parameters  Rcore  and  Rscaie •  The  values  for  the  parameter  Rcore  spanned  the 
range  0.02  -  1.0  kpc  and  the  values  for  the  parameter  Rscaie  spanned  the  range  0.005 
-  2.0  pc. 

In  addition  to  the  diffuse  models  described  above,  the  possibility  of  a  point  source 
of  511  keV  line  emission  in  the  direction  of  the  Galactic  center  was  also  investigated 
by  modeling  the  response  to  individual  point  sources  on  a  0.1°  X  0.1°  grid  within  the 
range  —5°  <  l  <  5°  and  —5°  <  b  <  5°.  A  point  source  located  at  the  position  of  the 
X-ray  source  IE  1740.7-2942  (359.14°,  —0.09°)  was  also  investigated  separately. 


ANALYSIS  AND  RESULTS 

The  surface-brightness  models  described  above  were  used  to  fit  the  narrow  511 
keV  line  data.  The  OSSE  data  fitted  include  those  reported  by  Purcell  et  al.  (1993) 
as  well  as  additional  data  collected  during  the  period  1992  November  17  -  1993 
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Figure  1:  The  OSSE  fields-of-view  for  the  observations  used  in  the  current  analysis, 
shown  in  Galactic  coordinates.  The  contours  represent  the  50%  response  levels  for 
each  of  the  OSSE  observations.  The  fields  used  as  background  measurements  are 
indicated  by  the  dashed  contours. 


September  7.  A  paper  describing  the  latter  data  is  in  preparation  (Purcell  et  al. 
1994).  The  OSSE  fields-of-view  for  the  observations  used  in  this  analysis  are  shown 
in  Figure  1.  The  OSSE  observations  of  the  Galactic  center  region  show  no  evidence 
for  significant  time- variability  of  the  511  keV  line  flux  (see  Figure  2  and  the  discussion 
below),  so  only  time-independent  models  were  investigated.  Including  the  results  from 
several  Gamma- Ray  Imaging  Spectrometer  (GRIS)  observations  of  the  Galactic  plane 
and  Galactic  center  were  also  investigated.  These  additional  data  included  the  results 
of  the  GRIS  observations  of  the  Galactic  center  in  1992  April/May  (Leventhal  et  al. 
1993),  since  they  were  nearly  coincident  with  one  of  the  OSSE  observations  of  the 
Galactic  center,  and  the  GRIS  observation  of  the  Galactic  plane  (Z  =  335°)  in  1988 
October  (Gehrels  et  al.  1991).  The  latter  was  investigated  since,  because  there  was 
little  response  to  the  Galactic  center  region  during  this  observation,  it  represents  a 
measure  of  the  diffuse  component  of  the  emission  independent  of  any  potential  point 
source  near  the  Galactic  center. 

For  each  separate  observation,  the  response  for  each  normalized  surface- 
brightness  model  was  calculated  by  convolving  the  2-dimensional  angular  distribution 
with  the  corresponding  instrument  angular  response.  The  OSSE  photopeak  angular 
response  used  for  this  analysis  was  based  on  detailed  angular  calibrations  performed 
prior  to  launch  (Johnson  et  al.  1993).  The  GRIS  angular  response  was  taken  to  be 
(Skibo  et  al.  1992): 


r(d) 


1  -  0.0570  0  <  16.75° 

0.067  -  0.00130  0  >  16.75° 


(7) 


For  the  OSSE  observations,  the  determination  of  the  model  response  included  the 
instrument  response  for  the  offset-pointed  background  positions  (see  the  discussion 
in  Purcell  et  al.  1993).  For  the  GRIS  observations,  it  was  assumed  that  the  back¬ 
ground  observations  were  taken  sufficiently  far  from  the  Galactic  plane  that  there 


Table  1:  Free  Spheroidal  Model  Fit  Results 


Dataset 

Model  Fluxa 

X2  Ndof  P(x2,  Ndof)  Spheroidal6  Disk6 

TOTALc 

OSSE 

OSSE  +  GRIS 

135.97  95  0.004  1.5  ±  1.0  0.5  ±  0.1 

136.72  98  0.006  1.7  ±0.6  0.5  ±  0.1 

2.0  ±0.9 
2.2  ±0.7 

a  (xlO-3  7  cm-2  s-1)  within  0.5  radian  of  l  =  0°,  b  =  0°. 

6  Uncertainties  represent  1<t  confidence  range  for  two  parameters  (Xmin  +  2.3). 
c  Uncertainties  represent  1<t  confidence  range  for  one  parameter  (Xmin  +  1-0). 


was  an  insignificant  contribution  from  the  source  distribution  during  the  background 
measurements. 

The  model  fluxes  were  determined  by  fitting  the  model  responses  to  the  data 
using  weighted  linear  least-squares  fitting  techniques  (Eadie  et  al.  1971).  The  covari¬ 
ance  matrix  used  in  the  fit  was  generated  using  the  statistical  uncertainties  in  the 
fitted  511  keV  line  fluxes  and  included  the  appropriate  correlation  terms  between  the 
OSSE  scanning  observations.  These  correlations  were  introduced  in  the  background 
estimation  process  where,  for  a  given  observation  period,  the  same  background  mea¬ 
surements  were  used  by  multiple  source  positions.  The  values  of  the  correlation 
coefficients  typically  ranged  between  —0.08  and  0.27,  with  a  median  value  of  ~  0.11. 

With  the  exception  of  the  Free  Spheroidal  Model,  all  of  the  models  tested  can 
be  rejected  at  a  high  level  of  significance  (^  4.5(r).  The  fitted  ffuxes  for  the  Free 
Spheroidal  Model  components  are  given  in  Table  1.  Of  the  models  investigated,  this 
model  provides  the  most  acceptable  fit  to  the  data.  The  M31  Novae  Model  provides 
the  next-best  fit  (%2  =  173.88,  P(%2,  98)  =  4  X  10— 6 ) ,  with  all  of  the  remaining  models 
having  probabilities  <  10-8.  Including  the  GRIS  observations  in  the  model  fits  does 
not  reject  any  additional  models,  but  does  provide  constraints  on  the  Free  Spheroidal 
Model  parameters  Rcore  and  Rscaie •  The  GRIS  results  were  included  in  the  model  fits 
discussed  below. 

There  is  no  evidence  for  significant  variability  of  the  511  keV  line  flux  during 
observations  of  the  Galactic  center  (see  Figure  2);  small  differences  in  the  observed 
ffuxes  can  be  explained  by  differences  in  the  observation  configuration.  For  the  Free 
Spheroidal  Model,  the  best-fit  values  for  Rcore  and  Rscaie  are  0.3  kpc  and  1.2  pc, 
respectively.  The  lcr  confidence  constraint  for  the  core  size  is  Rcore  <  0.5  kpc.  The 
Free  Spheroidal  Model  fit  is  not  improved  significantly  if  the  disk  component  given 
by  Equation  6a  is  replaced  by  any  of  the  other  disk  models  investigated,  indicating 
that  the  disk  component  of  the  emission  is  not  significantly  constrained  by  the  current 
observations.  The  range  of  fitted  ffuxes  for  different  disk  components  is  (3-5)  X  10“4 
7  cm-2  s_1  within  0.5  radian  of  l  =  0°,  b  =  0°. 

For  the  Point  Source  Grid  Model,  the  best-fit  point  source  position  is  found  to 
be  l  =  0.7°  ±  0.6°,  b  =  —0.8°  ±  0.3°,  and  is  consistent  with  the  Galactic  center  at  the 
~  3cr  level.  The  slight  offset  in  the  best-fit  point  source  position  provides  marginal 
evidence  for  an  asymmetry  in  the  distribution  of  the  511  keV  line  emission;  however, 
the  asymmetry  can  not  be  explained  by  assuming  additional  511  keV  line  emission 
from  the  x-ray  source  IE  1740.7-2942.  Adding  a  point  source  at  the  location  of 
IE  1740.7-2942  to  the  Free  Spheroidal  Model  does  not  significantly  improve  the  fit. 
The  95%  confidence  upper  limit  to  time-averaged  narrow  511  keV  line  emission  from 
IE  1740.7-2942  for  this  model  is  6  X  10-5  7  cm-2  s_1.  This  model-dependent  limit 


sion,  which  represents  ~  40%  —  60%  of  the  estimated  disk  component,  should  roughly 
follow  that  of  the  observed  26 A1  (it  is  interesting  to  note  that  the  COMPTEL  26 A1 
observations  suggest  the  central  peak  is  located  at  a  longitude  of  ~  2°  (Diehl  et  al. 
1994),  which  might  explain  the  slight  asymmetry  in  the  511  keV  line  distribution). 


Figure  2:  Averaged  narrow  511  keV  line  flux  from  the  OSSE  observations  within  ~  1° 
of  the  Galactic  center.  The  squares  indicate  the  expected  fluxes  for  the  best-fit  Free 
Spheroidal  Model. 
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Figure  3:  The  history  of  the  511  keV  line  flux  from  the  Galactic  center  region  as 
observed  by  various  balloon  and  satellite  instruments.  Note  that  the  uncertainties 
represent  the  statistical  uncertainties  only.  The  diffuse  contribution  for  the  best-fit 
Free  Spheroidal  Model  has  been  subtracted  from  the  reported  flux  for  each  observa¬ 
tion.  The  model  fit  included  all  of  the  OSSE  data  and  only  the  GRIS  data  points 
indicated  by  the  filled  circles.  The  numbers  in  parenthesis  indicate  the  field- of- view  of 
the  instrument,  in  degrees.  References:  SMM  -  (Share  et  al.  1990);  GRIS  -  (Gehrels 
et  al.  1991;  Leventhal  et  al.  1993);  Bell/Sandia  -  (Leventhal  et  al.  1978,  1980,  1982, 
1986);  FIGARO  -  (Neil  et  al.  1990);  HEXAGONE  -  (Chapuis  et  al.  1991);  CESR  - 
(Albernhe  et  al.  1981);  UNH  -  (Gardner  et  al.  1982);  GSFC  -  (Paciesas  et  al.  1982). 


The  remaining  positrons  are  probably  created  by  the  decay  of  56Co  and  44Sc,  which 
are  thought  to  be  produced  primarily  in  Type  la  supernovae  (Chan  and  Lingenfelter 
1993).  The  511  keV  line  flux  for  the  spheroid  component,  ~  1.7  X  10-3  7  cm-2  s_1, 
corresponds  to  a  positron  production  rate  of  ~  7  X  1042  positrons  s_1  assuming  the 
positrons  are  in  equilibrium.  If  the  positrons  are  primarily  produced  by  the  decay 
of  56Co  and  44Sc,  this  positron  production  rate  implies  a  limit  to  the  present  rate 
of  Galactic  iron  nucleosynthesis  of  ~  0.7  M@  per  100  years  (see  Ramaty,  Skibo  & 
Lingenfelter  1994  and  references  therein).  A  better  determination  of  the  distribution 
of  the  511  keV  line  emission  will  provide  more  significant  constraints  on  the  positron 
production  rate  for  the  disk  and  spheroidal  components,  and  on  the  present  rate  of 
Galactic  nucleosynthesis  of  /3+-unstable  nuclei. 
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